Complexity and heights of tori 

Gautam Chinta Jay Jorgenson Anders Karlsson * 
November 1, 2011 



Abstract 

We prove detailed asymptotics for the number of spanning trees, 
called complexity, for a general class of discrete tori as the parame- 
ters tend to infinity. The proof uses in particular certain ideas and 
techniques from an earlier paper [CJK10 . Our asymptotic formula 
provides a link between the complexity of these graphs and the height 
of associated real tori, and allows us to deduce some corollaries on the 
complexity thanks to certain results from analytic number theory. In 
this way we obtain a conjectural relationship between complexity and 
regular sphere packings. 

1 Introduction 

The number of spanning trees t(G), called the complexity, of a finite graph 
G is an invariant which is of interest in several sciences: network theory, sta- 
tistical physics, theoretical chemistry, etc. Via the well-known matrix-tree 
theorem of Kirchoff, the complexity equals the determinant of the combina- 
torial Laplacian Aq divided by the number of vertices. 

For compact Riemannian manifolds M there is an analogous invariant 
h(M), the height, defined as the negative of the logarithm of the zeta- 
regularized determinant of the Laplace-Beltrami operator, and which is of 
interest for quantum physics. The analogy between the height and complex- 
ity has been commented on by Sarnak in [S90J. 

In statistical physics it is of interest to study the asymptotics of the 
complexity, and other spectral invariants, for certain families of graphs. Im- 
portant cases to study are various subgraphs of the standard lattice Z d . An 
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instance of this is to study discrete tori, corresponding to periodic bound- 
ary conditions, as the parameters tend to infinity, see |DD88] . |CJK10| . and 
references therein. It is shown in |CJK10| that in the asymptotics of the 
complexity of discrete tori, the height of an associated real torus appears as 
a constant. 

In the present paper we study discrete tori of a more general type, defined 
as follows. Let A be an invertible r x r matrix with all entries being integers. 
This matrix defines a lattice AZ r in W. We associate to the group quotient 
with standard generators 

AZ r \Z r 

its Cayley graph, which we call a discrete torus. In other words, two elements 
x and y in AZ r \Z r are adjacent, denoted x ~ y, if they differ by ±1 in exactly 
one of the coordinates and equal everywhere else (everything mod AZ r of 
course). 

Let = Ao < Ai < ... < A | dot A|-i be the eigenvalues of the combinatorial 
Laplacian A^r^r of the discrete torus - see ([1]) for a definition of A^ r \Z r - 
Define det* A^r^f to be the product of the nonzero eigenvalues of the 
Laplacian: 

det*A AZ r\ Z r := AiA 2 ...A| de tA|-i- 

Note that the trivial eigenvalue is removed. We will nevertheless refer to 
det*A A z'-\Z'' as the determinant of the Laplacian. For simplicity we will 
mostly assume that det A > 0. 

Theorem 1 Let {A n } be a sequence of r x r integer matrices. Suppose that 
det A n —7- oo and A n /(det A n ) 1//r — > A 6 SL r (K). Then as n — > oo, 

2 

log det* A AnZ r\ Z r = c r det A n + - log det A n + log det* A AZ r\ Rr + o(l) 
where 

r°° dt 

c r = log2r- / e- 2rt (I (2t) r - 1)-. 

The definition of log de^A^r^r will be recalled in section 3 below. 

Our earlier paper |CJK10| treats the case when the A n are diagonal 
matrices. The present paper uses several facts that are established in that 
paper. Thanks to the fact that the first two terms in the asymptotics are 
universal in the sense that they only depend on A n via det A n , the theorem 
gives a close connection between the complexity of certain graphs and the 
height of an associated manifold. We emphasize that this is not at all obvious: 
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while it is true that the appropriately rescaled eigenvalues of the discrete 
tori converge to the eigenvalues of A^ r \]R r > this convergence is certainly 
not uniform. Moreover, the height cannot be defined as the product of 
eigenvalues, there is a regularization in the definition. An explicit expression 
for the height of flat tori of volume 1, h(AJ J r \M r ) := — log det A^rx Rr , can 
be found in Theorem 2.3 of [Ch97j. Dcningcr and Luck have informed us 
that the constant c r also has an interpretation as a determinant, namely the 
Fuglede-Kadison determinant of the Laplacian on Z r , see [L]. 

We turn now to a connection between our results and sphere packings. 
The problem of finding the densest packing of ordinary space with spheres 
of equal radii is an old one with practical importance even in dimensions 
greater than 3. One type of packings is regular sphere packings which means 
that the spheres are centered at the points of a lattice AT? . Gauss showed 
that the face centered lattice (fee) D% is optimal among regular packings in 
three dimensions. By work of Thue and Toth one knows that the hexagonal 
lattice Ai is densest in dimension 2. In dimenison 24 it is known that the 
Leech lattice is optimal among unimodular lattices. We refer to the book of 
Conway-Sloane [CS99J for more information. 

Conjecturally the height of ^4Z r \IR r has a global minimum when All is 
the densest regular sphere packing. Extremal metrics for heights has been 
studied in [OPS88J in dimension 2, and for tori in higher dimensions notably 
in |Ch97j and |SS06] , In these papers, the question is phrased as the study 
of the derivative of Epstein zeta functions at s = 0. From this theory we can 
deduce the following corollaries from our main theorem: 

Corollary 2 Let {A n } be a sequence ofrxr integer matrices with det A n — > 
oo. Suppose that {A n /(det A n ) 1 / r } belongs to a compact subset of SL r (M), 
r = 2,3, avoiding lattices equivalent to A2, resp. D3. Assume that there is 
a sequence {L n } with detL n = det A n such that {L n /(det L n ) 1 / r } converges 
to A2, resp. D3. Then L n 7?\T? has more spanning trees than A n r E r \L r for 
all sufficiently large n. 

Corollary 3 Let A n be a sequence ofrxr integer matrices with det A n — > 00 
Suppose that {A n /(det A n ) 1//r } stays in a compact subset of SL r (M). For all 
sufficiently large n we have that 

r(A n Z r \Z r ) < 1 " j exp(c r det A n + 7 + 2/r), 

47T 

where 7 is Euler's constant and c r is as in the theorem. 
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In the trivial case r = 1, this estimate gives a value close to the truth: 



det A n = T{A n Z r \Z r ) < 1.05 det A. 



Upper bounds for the number of spanning trees have been considered in 
the combinatorics literature since 1970s at least. For regular graphs there is 
a rather sharp estimate by Chung and Yau |CY99| improving on an earlier 
result of McKay [Mj. In general, it is an open problem to decide which 
simple graph on n vertices and e edges has the maximal complexity. This 
is of interest to communication network theory since this graph invariant 
appears as a measure of reliability. 

It would be of interest to also go in the other direction: proving results on 
the extrema of families of Epstein zeta functions via a better understanding 
of the number of spanning trees of discrete tori. 

2 Spectral preliminaries for discrete tori 

Let A be an invertible r x r integer matrix. This matrix defines a lattice 
AZ r in W . We denote by DT(A) the discrete torus, or Cayley graph of the 
quotient group AZ r \27 with standard generating set: two elements x and y 
in AZ r \Z r are adjacent, denoted x ~ y, if they differ by ±1 in exactly one of 
the coordinates and equal everywhere else (everything mod AZ r of course). 
The associated (combinatorial) Laplacian is defined by 



on functions / : AZ r \Z r -)• R. 

The dual lattice A*Z r is as usual all the points v in W such that (v, x) £ Z 
for all x £ AW, where (•,•) denotes the usual scalar product. Since W is 
self-dual and AW is a subgroup of Z r it follows that 27 is a subgroup of 
A* 27. Note that the respective indices are 





(1) 



y s.t. i/~a; 



[IT : AW] = [A* :Z7] = |detA|. 




r 
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where vt denotes the kth coordinate of v. 

Proof. The operator AjxrfA) i s a semipositive, symmetric matrix and hence 
we are looking for |detA| number of eigenfunctions and eigenvalues. The 
proof is a trivial calculation: 



) 



z~0 



The heat kernel K A (t,x) : R>o x DT(A) — >• R is the unique bounded 
function which satisfies 



A DT (A) + ^) K A (t,x) = 

K A (0,x) = 5 (x), 



where 8o{x) = 1 if x = and otherwise. The existence and uniqueness of 
heat kernels in a general graph setting is established in [DM06] . Recall from 
e.g. |C.TK10j that 

r 

K zr (t,z) = l[K z (t,z k ) 

k=l 

for z = (zfc) and K z (t,w) = e~ 2t I w (2t), where I w is the 7-Bessel function of 
order w. 

We have the following theta inversion formula (cf. |CJK10| ). 
Proposition 5 The following formula holds for x £ AZ r \Z r and t £ R>o 

j/eAZ r 1 1 vez r \A*z r 

In particular, 

9 A (t):=\detA\ £ e - 2ri /, 1 (2t).../, r (2t) = £ 

Proof. Since both sides of the equation satisfy the conditions for being the 
heat kernel, this follows from the uniqueness of heat kernels. The second 
formula is the special case x = 0. ■ 
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3 Spectral preliminaries for continuous tori 

An r- dimensional (continuous) torus is given as a quotient of W by a lattice 
AH where A 6 GL r (M). The metric structure and the standard (positive) 
Laplace-Beltrami operator — *}2id 2 /dx 2 on ^ r projects to the torus. The 
volume is |det A\ . Let A* be the matrix defining the dual lattice A*17 , and so 
A* = (^4 _1 ) • The eigenfunctions of the Laplacian on the torus in question 
are f v {x) = exp(2iriv t x) where v are the vectors in the dual lattice. The 

9 1 1 1 1 2 

corresponding eigenvalues are = 4ir \\v\\ or with a different indexing: 
X m = (2ir) 2 (A*m) t (A*m), where m runs through 17 . We have the associated 
theta function 



meZ r 

The theta inversion formula, which is equivalent to the Poisson summation 
formula in this case, yields 



The associated spectral zeta function, which in this case also goes under 
the name of the Epstein zeta function, is defined as 



Classically, one can prove the meromorphic continuation of Z\{s) to all 
s G C, showing that its continuation is holomorphic at s = 0. From this, one 
defines the spectral determinant det*A j 4 Z ryj R r by 






log det*A AZ r\ R r 



Z' A (0). 



4 Asymptotics 



Let 




and 



Ha{s) = - / ( e~ sH [9 A (t) - |det A| • e~ 2rt I ( 



o 



(2tY - 1] + e"«) ^ 



Jo 
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Everything in section 3 of [CJKlOj carries over with only notational changes, 
even though the eigenvalues are different and the theta identity is hence 
somewhat different. These differences are not essentially used. In particular 
the first order term as t — > in the trace of the heat kernel is still (in 
the present notation) | det A| • e~ 2rt Io(2t) r since it corresponds to the trivial 
eigenvalue. In particular the following extension of Theorem 3.6 in [CJKlOj 
holds: 

Theorem 6 For any s £ C with Re(s 2 ) > 0, we have the relation 
log(s 2 + A^) = |det A| -X r (s) +H A (s). 

Letting s — > we have the identity 

log( Yl A,) = | det A| -2^(0) +Ha{0). 

Section 4 of |CJK10| is an independent section on uniform bounds on 
/-Bessel functions. We recall the following statements, slightly adapted to 
the present context (keeping in mind that I- y = I y for integers, and that b 
may now be real): 

Proposition 7 For any t > and b > 7 there is a constant C such that 

< V¥te- h2t h{b 2 t) < C < 1 

Fix t > and integers y, no > 0. Then for all b > uq we have the uniform 
bound 

-n y/2b 



0<VfcW^ t )<( 1 + _y 



Proposition 8 Let N(u) be a sequence of positive integers parametrized by 
u £ Z + such that N(u)/u — > a > as u — > oo. Then for any t > and 
integer k, we have 

hm N(u)e^ h L N{u)k {2uh) = " e~^^. 

From now on we fix a sequence {A n } of integer matrices with < 
det A n — > oo satisfying 

■ — — -rA n — > A as n — > oo, 

for some A £ SL r (R). From the previous propositions we will deduce the 
following: 
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Proposition 9 For each fixed t > we have the pointwise convergence 

A „(det(A n ) 2 / r i) -> 6 A {t) 

as n —> oo. 

Proof. For any »eZ r and A G GL r (R), let 



i=l 



where (At>)j denotes the i-th component of Aw. Let u n = det(A n ) 1 / r and 
cii = (Av) i . Note that (Av)i/u n —> <2j. We have 



From Proposition [8] (with fc = or ±1) we have for any t > and u 6 
that 



as n — )• oo. This means that the proposition will be proved if we can inter- 
change the limit and the infinite sum. We show that for fixed t, the sum is 
convergent uniformly in u n (or equivalently, in n). 

We can rewrite the sum 9\ n (u^t) in r + 1 sums depending on how many 
components of the A n v are zero. Pick no sufficiently large so that 

K| /2 < \(A n v)i\ /u n < 2 \ai\ 

for all !)£Z r and n> uq and dj 7^ 0. Recall that I_ n = I n . Let us look at 
a term with k zeros in the yis and estimate with the help of Proposition 



/ -1 \r / |/a \ i\ -no\(A n v)i\/2u n 

U) r n Akl 



r 

< 1 1 TT \l a »l 



for all n large and where 

\-no/4 
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and a is the smallest nonzero absolute value of all the entries in All . The 
whole theta series is therefore bounded by r + 1 sums of a product of con- 
vergent geometric series. This shows that the infinite sum is uniformly con- 
vergent and the proof is complete. ■ 

Lemma 10 Given a sequence {A n } satisfying (det A n ) _1 / r A n — > A as above, 
there is a constant d > such that for all sufficiently large n 

oo 

*A n (t#)<l + jV«*. 

i=i 

Proof. Let u n = det(A n ) 1 / r . We have 

9 An (t)= e ~ tXv = Yl e"*( 2r " 2E ^i cos(27 ™ fc) ) 

r 

= 1 + Y Yle~ 4ts ' m2{wVk) 

v&L r \h*JL r k=l 

So that 

r 

fl An (n 2 t)=i+ n e ~ 4tu " sin2( ™ fc) 

ueZ r \A^Z r k=l 

We use the elementary bounds sinx > x — x 3 /6 and sin(7r — x) > x — x 3 /Q 
for x G [0,7r/2] and get 

u n sm(nvk) > u n irv k (l - 7r 2 ^/6) > ciru n v k if v k < 1/2 
■u n sin(7pu fc ) > u n irv k (l - Tr 2 v k /6) > ciru n (l - v k ) if v k > 1/2 

for some postitve constant c for all n sufficiently large. Note also that for 
every v ^ 0, the values u n v k range over the integers times an entry in A as 
n — > oo because of the convergence of (det A n ) -1 / r A n to A in SL r (R). We 
then conclude there is a constant d > such that for all sufficiently large n 

oo 

An (t#)<l+I>" dy - 

■ 

Now we can show: 
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Proposition 11 With the notation as above and u n := det(A n ) 1//r , we have 
that 



/DC . 
(e An (u 2 n t) - ule- 2 ^h{2ulty -1 + e- 

(e A (t)-i) T --(^r'/ 2 + o(i) 



t\ dt 

t 



as n — > oo. 



Proof. Write 



(0 A >*i) - <e- 2 ™«%(2^r - 1 + e-< 1 
= f (0A n «t) -1) j- f u r n e-^I (2ultY^ + jf' 



-u 2 1 ^ 

e " — . 
t 



In the last row, the third integral clearly goes to zero as n — > oo. For the 
first integral in the same row we have 

in view of the pointwise convergence from Proposition [9] and the uniform 
integrable upper bound from Lemma [TU1 
The middle integral 



/oo 



2 j\r dt 



converges to 

\ 47Tt) -r/2dt = ^ (4vr) -r/2 

in view of the heat kernel convergence from Proposition [7] and Proposition 
El so then we may appeal to the Lebesgue dominated convergence theorem. 
■ 

Next we show: 

Proposition 12 With the notation as above and u n := det(A n ) 1//r , we have 
that 

dt f /„ , ^ ri ^_wo\ dt 

T 



(e An (u 2 n t) - uie- 2 ^ t i {2ulty) a - -+ / (e A (t) - ^r 72 



as n — > oo. 
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Proof. For fixed t we have the pointwise convergence as n — > oo 

9 An (u 2 n t) - <e- 2 ™«%(2^r -> A (*) - (4vrt)- r / 2 . 

It remains therefore to exhibit uniform (for n » 1) integrable bounds on the 
integrands. This can be done in the same way as in the proof of Proposition 
[9j In order to make the bound obtained there, in terms of A, integrable for 
< t < 1, we just need to choose no large so that noa/4 > r/2. ■ 
Finally we recall from |CJK10| : 

Proposition 13 For u £ M. we have the asymptotic formula 

/' 1 ( e - 2 t_l ) ^ = r'(l)-21og(u)+o(l) 
Jo T 

as u —> oo. 



We now turn to the proof of our main result, Theorem 1. 
In view of Theorem [6] we have 

/•oo 

logdet A DT{An) = det K-c r - / {e An (t) - det A n • e - 2rt I {2t) r - 1 + e"*) 

J o 

After the change of variables t — > u^t, the second term becomes 
- (0A n (u 2 n t) - det A n • e-^Wulty - 1 + e"^*) ^ 

f°°l / 9 9 \ fit 

+ / ( 0A n «t) ~ det A n • e- 2 ™«*I (2^i) f - - 1 + e' u A -. 



dt 
T' 



In view of Propositions QTJ (T2J and [T3l this integral equals 



°°(e A (t)-i)^ + k^r r/2 

t r 



r'(l)+21og(n n ) + (l). 



1 A (t) - (47Ti)- ? -/ 2 ) 



dt 

T 



Keeping in mind that u n = (detA n ) 1//r and identifying the constant terms 
appearing in the meromorphic continuation of —£'(()) = logdet A-AZ r \R r i the 
main theorem is proved; cf. equation (15) of |CJK10] with V(A) = det A = 1. 
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5 Proof of the corollaries 



To prove the corollaries in the introduction we recall the statements from 
the literature that we use. 

For Corollary [2] note that the height has a global minimum for the hexag- 
onal lattice A2 in dimension 2 as is well-known and for the f.c.c. lattice 
A3 = D3 in dimension 3 by the rigorous numerics of Sarnak-Strombergsson 
in |SS06| . Hence for any unimodular lattice L in the respective dimensions 

h(L) > h(A 2 ), so logdet A L < logdet A Ai 
h(L) > h(D 3 ), so logdet A L < logdet A D:i . 

As already remarked the two leading terms in the asymptotics in Theorem [1] 
are shape-independent, and so Corollary[2]follows from these remarks arguing 
with convergent subsequences in view of the compactness. 

For Corollary [3] note that Corollary 1 on p. 119 in |SS06] implies that 

2 

log det Am < 7 - log 47r + - < -0.95 

r 

where M is an r-dimensional flat torus of volume 1 and 7 is Euler's constant 
7 ~ 0.577. In view of this statement, let us replace A n by a convergent 
subsequence. Then by Theorem Q] and the matrix-tree theorem, we have 

r(DT(A n) ) = ^ < (det K)^ . + + 

det A n 47r 

for all sufficiently large n. This concludes the proof of Corollary [3l 

Finally, it may be of interest to mention another estimate in [SS06J: 

^vK"§f) r(i+o(i)) 

where m(L) is the length of the shortest non-zero vector in the lattice L. 
Recall that being the densest regular packing is equivalent to being the lattice 
with co- volume 1 which maximizes the length of the shortest nonzero vector. 
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